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In the CMR materials, ferromagnetic double exchange (DE) presumably coexists with a direct nearest- 
neighbour antiferromagnetic interaction. We construct a single-site mean field theory that explicitly takes 
into account the different nature of carrier-mediated ferromagnetism vs. Heisenberg-like superexchange. 
We find, in contrast to previous results in the literature, that the competition between these two exchange 
interactions leads to ferro- or antiferromagnetic order with incomplete saturation of the magnetization (or 
sublattice magnetization), rather than spin canting. The associated experimental implications are discussed. 
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I. INTRODUCTION 

Recently, there has been a renewed interest (moti- 
vated by technological problems of microelectronics) in 
the propertieSpOf colossal magnetoresistance (CMR) man- 
ganese oxidesEl . The CMR behaviour typically corre- 
sponds to an intermediate doping range, when these ma- 
terials are ferromagnetic. The latter property is gener- 
ally attributed to a conduction electron-mediated dou- 
ble exchapge (DE) interaction^. In addition, there exists 
evidencetTEl which suggests the presence of antiferromag- 
netic superexchange interactions of comparable scale. In 
this paper, we study the behaviour of a classical magnet 
with competing double exchange and superexchange in- 
teractions, and show, in particular, that in an isotropic 
case the spin canting (which was previously suggestedtl 
to be a generic outcome of such a competition) can be 
stabilized only at very high fields and at very low tem- 
peratures. |-| 

Previous related calculationsO have been performed 
for a strongly anisotropic model in which the inter- 
and intrajlayer direct exchange constants have differ- 
ent signsfl. Here, we assume that direct interactions 
have everywhere the same (antiferromagnetic) sign and 
magnitude. This is viewed as more appropriate for the 
Lai-xCaxMnOs perovskite family away from the a; = 
endpoint, as well as for the layered manganates such 
as La2-2a;Sri+2a;Mn207. In these compounds, pairs of 
stacked Mn-0 planes form the bilayers, which are sepa- 
rated by poorly conducting layers of La(Sr)0 (see Ref. |^). 
Since the lengths of the intra- and inter-layer bonds in 
a bilayer are roughly the same, the values of inter-layer 
hopping coefficient and superexchange constant should 
be of the same order of magnitude as their intra-layer 



counterparty. The observed interlayer cantingEI or canted 
correlationso, presumably caused by superexchange be- 
tween the two layers of a bilayer complex, can be used 
to suggest relatively large values of superexchange within 
the layers as well. In view of the considerable interest in 
these layered systems, we consider primarily the two di- 
mensional (2D) lattice. By doing so, we expect to capture 
the basic magnetic properties of the layered compounds, 
while avoiding the cumbersome quantitative treatment 
of bilayers. We also note that our 2D results are qual- 
itatively representative of the three dimensional case as 
well. 

We will see that the fact that the double exchange ~ 
superexchange competition occurs at all lattice bonds (i. 
e., in the 3D case, for in-, as well as out-of-plane inter- 
actions) leads to the enhancement of spin fluctuations. 
Lowest order (i.e., Hartree-Fock like) treatments arc in- 
sufficient in this case. Therefore, in the present paper we 
introduce a new approach to the problem. 

We construct a single-site mean field theory which ex- 
plicitly takes into account the main feature of the present 
problem, namely, the carrier kinetic energy origin of the 
double exchange ferromagnetism. We note that thejaean 
field schemes previously reported in the literatureQiJ es- 
sentially use an effective Heisenberg-like ferromagnetic 
exchange interaction to describe the double exchange. 
Such schemes would not adequately reflect the very dif- 
ferent nature of the two competing interactions. It is 
of interest, then, to see how the results are changed if 
a more proper treatment of the double exchange is |Cai:- 
ried out. We begin with the standard HamiltonianHEj, 
corresponding to an infinite on-site Hund's rule coupling: 
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Here, the first term is the kinetic energy of the carri- 
ers (which are represented by the fermion operators Cj 
where j is the site index). The second term corresponds 
to the nearest-neighbour antiferromagnetic {Jaf > 0) 
exchange interaction between the classical {S ^ 1) core 
(localized) spins Si, and the last term is the Zeeman en- 
ergy of these spins in the external field H. The double 
exchange interaction results in the modulation of carrier 
hopping coefficients by the factors cos0.y = Si ■ Sj/S'^. 
Since in the present work we restrict ourselves to a single- 
site mean field treatment, we omit the additional phase 
factors which would result in Berry phase effects. We use 
units in which the bare hopping coefRcient to, h, ks, fJ-B, 
and the lattice spacing are all equal to unity. 

Our mean field framework is based on the exact solu- 
tion of the single site problem, which is outlined in the 
following section. Details of the derivation are relegated 
to Appendix whereas the implications for the T = 
energetics are briefly discussed in Appendix ^ Section 
III is concerned with the mean field treatment of dif- 
ferent magnetic phases of the system, and the resultant 
mean field phase diagram is described in Section FV . We 
conclude with a brief |discussion of the experimental rel- 
evance of our findingsEJ. 



II. EXACT SOLUTION OF THE SINGLE-SITE 
PROBLEM 

The random distribution of localized spins leads in 
Eqn. (^ to a highly disordered electronic hopping prob- 
lem. Our mean field treatment is based on evaluating 
the energy cost, SF, of a fluctuation of a single spin 5*1, 
embedded in an effective medium with a uniform aver- 
age value of cos 6'y . As we will see below (Eqns. (^),(13), 
and ([l9|) ), such a fluctuation gives rise to a difference 



between the values of hopping coefficient b from the site 
of the fluctuating spin 5*1 to the neighbouring sites and 
the backgrpund hopping t {t b); for clarity these 
parametersla are indicated schematically in Fig. 1. 

The quantities b and t depend on in a self consistent 
fashion on the change, 6Foe, in the free energy, associ- 
ated with the local change in hopping matrix elements 
t ^ b. Such a local change, originating from a local spin 
fluctuation on the site (0, 0), gives rise to a perturbation, 
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i(6 - t) |cJq_„j (c(o,i) + C(i,o) + C(o,_i)-h 



+ C(_i.o)) +h.c.} 
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in the carrier kinetic energy. This perturbation shifts 
the energy levels of individual carriers, thus resulting in 
a change in the total kinetic energy of the band. This 
kinetic energy contribution to 5F , which can be evaluated 
following Refs. |l3|-|l5| (see Appendix ^) , is given by 



5FDE{b,t,T)^ / /(e)C(e)de 



+e{b-t)-{if{zo)~^{-Dt)) , 

where the spectral shift function ^(e) takes values be- 
tween — 1 and 1, and is given by 



C(e) = Arg jfo^ + (f - b')e j V 



v(ri)dr\ 
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the bound state energy zq < —Di is the root of 
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iy{e) is the density of states. 
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and fi is the chemical potential. Although we will apply 
Eqns. (||-H) only to the case of a 2D square lattice, they 
remain valid for a cubic lattice in 3D, as well as in the 
ID case. The energy integrations are performed over the 
conduction band width, —Dt < e,?/ < Dt where D is 
the dimensionality of the system. It should be stressed 
that it is because of the locality of the perturbation (||), 
which represents a lattice analogue of an s-wave scat- 
tering proUem, that the quantity SF can be evaluated 

The second term in Eqn. (|^) is a contribution of a 
bound state that appears in the carrier spectrum for b > t 
(when the perturbation may be viewed as a "potential 
well" ) . In 2IX-|the binding energy of this state vanishes 
exponentiallyO as 6 ^ i, whereas in 3D it has a threshold 
behaviour. This bound state is related, to one that causes 
the formation of magnetic polaronalS, but it should be 
distinguished from the true magnetic polaron which is an 
extended object and can not be treated within a single- 
site approach. 

In order togain additional intuition about the meaning 
of Eqns. (|3|-||), it is useful to calculate the energy cost of 
a single-spin fluctuation in various phases at T = 0. This 
is discussed in Appendix^. This Appendix highlights the 
important differences between the double exchange and 
familiar Heisenberg direct exchange interactions. 

For the purposes of the present work, the virtual crys- 
tal approximation, based on the parameters shown in Fig. 
1, is expected to be appropriate as long as the carrier 
concentration is not too small. This is because the quan- 
tities of interest involve integration over carrier energies 
in the metallic phase. In order to extend this formulation 
beyond the single site mean field scheme, we note that 
multi-site spin fluctuations can also be treated as local 
perturbations following Ref. ^ This in principle allows 
one to study systematically the effects of correlations, by 
constructing an analogue of an impurity-concentration 
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expansion. This procedure would also verify whether the 
virtual crystal approximation is a good starting point for 
studying other (e. g. transport) properties. 



III. THE MEAN FIELD SCHEME 

1. Ferromagnetic phase. In the ferromagnetic phase 
at r > 0, the net energy cost of a single-spin fluctuation 
is (in 2D) given by 

SFi = 6FDE{b,t,T) + 4:Jaf{cos9i2)2 - H"cosai - 

-4Jaf(cos6'i2)i2 +i?(cosai)i . (7) 

Here, 6112 is the angle between the directions of the fluctu- 
ating spin Si and any of The neighbouring spins, denoted 
by 5*2 (we assume that spin fluctuations on different sites 
are statistically independent), and ai is the angle be- 
tween Si and the direction of magnetization, M (see Fig. 
|l|). The angular brackets, (...)/, are used to denote the 
average values over the Boltzmann probability distribu- 
tion of spin Si, wi (X expL^SFi/T). We then flnd that 
(cos0i2)2 — Mcosai, andll3 



= (cos2(0i2/2))2 = (1 + Mcosai)/2, 
= {9)i = (l + Af2)/2. 



(8) 



The magnetization has to be determined self-consistently 
as M = (cosai)i; generally, the latter equation has to 
be solved numerically. 

In the ferro- and antiferromagnetic phases, it is use- 
ful to construct a reference framework with which to 
compare our results. We define Jeff(M) which repre- 
sents an effective M- dependent exchange constant for a 
Heisenberg-like magnet. The appropriate exchange con- 
stant can be deduced by considering small spin fluctua- 
tions (I cos a — A/ 1 <C 1), which correspond to small fluc- 
tuations in the hopping matrix elements ( |t — 6| <C i). A 
perturbation expansion of Eqn. (0) then leads to 
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ef{e)y{e)de = 2{t~b)\E^\, 



(9) 

at leading ordei@ in T/t, where i?o is the kinetic energy 
of the carriers for f = 1. In the ferromagnetic state, it 
follows from Eqn. (||) that 



1, 



Jtff{M) = JAF~-AE, 



1 + M2 



(10) 



The second term in the above equation represents the 
DE contribution. This-.term, which is contained in other 
mean field schemesoO, increases as M decreases. As 
a consequence, for moderately strong antiferromagnetic 
exchange interactions, whenE3 



J^^j (M) changes sign as M varies from to 1. This be- 
haviour has important consequences: it leads to a lack of 
saturation in the low temperature magnetization. Typi- 
cal results for M{T) are plotted in Fig. ^ for these mod- 
erately strong exchange interactions. Here the solid line 
represents the full mean field calculation (which makes 
use of Eqns. (§-§)), while the dashed line corresponds to 
the effective exchange approximation. The dotted line 
represents the behaviour of a conventional Heisenberg 
magnet with the same value of Curie temperature, and a 
constant nearest-neighbour exchange integral. 

The lack of saturation seen in Fig. || can be under- 
stood as follows. In the paramagnetic phase, M = 
and, by virtue of Eqns. (|Tc|-ll), the effective exchange 
constant has a negative (ferromagnetic) sign. As T 
decreases, the system undergoes a Curie transition at 
Tc ~ 4|Jg^y (0)1/3. Decreasing T further results in a 
decrease in the magnitude of spin fluctuations, i. e., in 
an increase of M. The latter is opposed by a decrease in 
I Jfjr^^l, leading to the softening of spin fluctuations. As 
a result, the effective exchange constant "self- adjusts" in 
such a way that it never becomes large in comparison 
with T, and even at low T the behaviour of an effec- 
tive exchange magnet is similar to that of a conventional 
Heisenberg magnet in a "high-temperature" regime of 
T Tc J. In this way, the non-vanishing thermal 
fluctuations do not allow the magnetization to reach its 
proper saturation value, Mq = 1. At zero flelcL4he value 
of magnetization as T ^ is instead given bynJ 



Mo = ^{Eo/JAFYIi2 - 1 < 1. 



(12) 



|So| < SJaf < V2\Eo\ , 



(11) 



These self consistent changes in | J^/j \ lead to inadequa- 
cies of the effective exchange approximation at low T. 
As may be seen in Fig. |^, the behaviour obtained in this 
approximation differs signiflcantly from that found using 
the full calculation of M (T) . This difference is due to the 
fact that when J^/f ^ T is small (in comparison with the 
electronic energy scales), quadratic terms (in {t — b)/t) 
dominate the physics. The details are discussed in Ap- 
pendix Within the effective exchange approximation, 
strong fluctuations of both angular co-ordinates of each 
spin persist at low T. By contrast, the full calculation 
shows that the fluctuations of the polar angle freeze out, 
cosai — > Afo- Independent fluctuations of the azimuthal 
angles of the spins, which persist down to T — > 0, appear 
to be an artefact of the single-site mean field treatment. 
It is natural to expect that, at least in the classical case 
of 5 ^ 1, these azimuthal fluctuations also freeze out (al- 
beit at a lower temperature than the polar ones) resulting 
in the formation of a multi-sublattice or spin-glass-like 
state with the net magnetization given approximately by 
Eqn. (^. 

2. Antiferromagnetic phase. The Neel antiferromag- 
netic state (of the metallic phase) can be treated sim- 
ilarly. We find (cos0i2)2 = —mcosai, where ai is the 
angle formed by the spin 5*1 with its average direction and 
m is sublattice magnetization. Eqns. (S) are replaced by 
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62 = (1 _ mcosai)/2, ^{\-m^)l2. (13) 
Instead of ([l^) we obtain 
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(14) 



It is easy to show that the Neel ordering arises for Jaf > 
2-5/2 |£;o I 2cro field and at T < Tat « U^f/'{0)/3. It 
always exhibits undersaturation of the sublattice magne- 
tization: at T — > 0, 
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= V1-(^^o/Jaf)732<1. 



(15) 



This undersaturation (which leads to a finite bandwidth) 
may be viewed as consistent with the presumed metallic 
state. 

Numerical calculations yield the dependence of m on 
T, which is similar to M{T) in the ferromagnetic phase 
and shows the same low-temperature features. 
3. Canted phase. Our discussion thus far has not io.- 
cluded the canted phase first proposed by De GennesQ. 
In our case, this is a two-sublattice (checkerboard) mag- 
netic phase; the sublattice magnetizations have an equal 
magnitude to and form an angle 27 with each other. In 
the present model, spin canting requires the presence of a 
magnetic field to break the high degeneracy which would 
otherwise occur. This degeneracy is related to the fact 
that the energy of the system depends solely on the values 
of the angles formed by the pairs of neighbouring spins. 
All the neighbours of any spin Si of sublattice I belong 
to sublattice II, and are parallel to each other at T = 0. 
Therefore, the energy of the system does not change as 
the spin Si moves along any cone around their common 
direction. In the context of single site mean field ap- 
proaches, the same holds at T > for any cone around 
the average direction of the sublattice II spins. Thus, 
the probability distribution of the spin will be axially- 
symmetric with respect to the direction of the magneti- 
zation of sublattice II, with which the spin Si will there- 
fore be aligned on average (rather than with sublattice I) . 
Thus, in the absence of perturbations (caused by next- 
nearest-neighbour exchange, anisotropy effects, quantum 
corrections, or small external fields) the canted statej-is 
destabilized, as a repilt of the underlying degeneracyo. 
Since it is site-loca£3, its effects will persist as long as 
the energy scale of a perturbation per individual spin re- 
mains small in comparison with the characteristic energy, 
/csT, of the thermal motion of a single spin. 

To characterize the finite field canted state, we use 
the full non-perturbative expression (|^). The mean field 
framework of Eqns. (0^^) has to be modified to allow 
for a self-consistent determination of the two mean-field 
variables, to and 7. We now obtain two coupled mean 
field equations, which, as in the ferromagnetic phase, fol- 
low from the self-consistent definition of the sublattice 
magnetization m. For the component of (5*1) parallel to 
the magnetization of sublattice I, we obtain 



— sin 27 (sinai cos Pi)i -{- cos 27 (cosq!i)i = to , (16) 
whereas the perpendicular component must vanish, 

cos 27 (sin «! cos/3i)i 4- sin27 (cosai)i =0. (17) 



In writing Eqns. (^-17) we used a co-ordinate system 
with a polar axis parallel to the sublattice II magnetiza- 
tion, cki and (3i are polar and azimuthal angles of the spin 
Si in this frame, with /3i = corresponding to the spin 
Si lying within the plane containing the two sublattice 
magnetizations. 

For the net energy of a single-site fluctuation we now 
obtain, instead of Eqn. (|^), 

5Fi= SFjjEibjtjT) + 4J^i?TOC0sai — 4:JAF'm{cosai)i — 
—H{— sin 7 sincki cos j3i + cos 7 cos ai) -I- 
+if(— sin7 (sinai cos/3i)i + cos7(cosai)i ) (18) 

whereas the values of the hopping coefficients (see Fig. 1 
and Eqn. (^) are given by 

6^ = (1 + mcosai)/2 , ^ (1 + cos27)/2 . (19) 

The low-T canted state is found to be stable for ^Jaf > 
\Eo\+H. 

We begin with the case of relatively large bandfilling, 
corresponding to the undersaturated ferromagnetic be- 
haviour at = (see Eqn. (|l|) and Fig. |). The 
solutiongSa of Eqns. (p^-[T7|) for typical parameters are 
illustrated in Fig. ||. One can see that, as T in the 
canted phase, the sublattice magnetization m approaches 
its proper saturation value m — 1. Note that the ferro- 
magnetic (7 = 0) solution to the mean field equations is 
present at > as well. In Fig. ^, the corresponding 
magnetization, Mfm{T), is represented by the dotted 
line. However, when the canted (7 > 0) solution exists, 
it corresponds to a lower value of the free energy. This is 
obvious from the fact that the net magnetization in the 
canted state Mcm{T) = —dF/dH = to cos 7 (dashed 
line in Fig. ^ is larger than Mfm{T)- The canted solu- 
tion branches from the ferromagnetic one at temperature 
Ti ^ H, when 



ATM = H{sm^ai)i 



(20) 



at this point the undersaturated ferromagnetic state un- 
dergoes a second-order spki-flip transition into the low 
temperature canted stateEII. One can therefore conclude 
that undersaturation is representative of the generic low- 
temperature b^iaviour of a double exchange - superex- 
change magneB^. 

For smaller values of carper concentration, at iJ > 
we find the spin-flop phaseEa of the undersaturated an- 
tiferromagnet, which evolves into the canted state via 
a smooth crossover at T ~ i?, as shown in Fig. ^ The 
low-temperature region where the canting angle 7 rapidly 
increases with T corresponds to the canted phase. 
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IV. MEAN FIELD PHASE DIAGRAMS 

Typical phase diagrams for the DE-superexchange 
magnet in (a,c) zero and (b,d) non-zero field are pre- 
sented in Fig. |5[ For io of the order of an eV, our 
choice of parameters corresponds to reasonable values of 
Jaf ~ 300K. In zero field (a,c), the solid line represents 
the phase boundary between paramagnetic (PM) and 
antiferro- (AFM) or ferromagnetic (FM) metallic phases. 
For the values of parameters used in Fig. || (a), the or- 
dered phases are undersaturated at low T. For slightly 
smaller Jaf we find a critical value of bandfilling, xi, 
which divides the saturated, x > xi, and undersaturated 
regimes (see Fig. ^ (c)). At low temperatures and small 
concentrations, the undersaturated AFM state becomes 
thermodynamically unstable {dfi/dx < 0), signalling ei- 
ther the onset of a more complicated spin arrangement 
or phase separation (see Appendix The dashed line 
in Fig. 1^ (a,c) corresponds to the anticipated boundary 
of this region {dfi/dx = 0). We note that the possibil- 
ity of phase separation in DE-superexchaiige systems has 
been suggestecL-hoth by analytical studies^ and numeri- 
cal simulations&ll. 

Figs. H (b,d) show that in the presence of a mag- 
netic field the PM-FM transition is replaced by a smooth 
crossover (dotted line). The spin arrangement of the 
AFM phase becomes non-coUinear (flop-phase), and has 
the same symmetry properties as the canted phase (CM), 
which becomes stable at lower T (replacing the H ^ Q 
undersaturated FM and AFM phases). The two are sep- 
arated from the PM and FM region by a second-order 
phase transition at T = Ti{x), which is represented 
by the solid line. At sufficiently small x the latter ap- 
proaches the H = Neel transition line. The thermody- 
namic instability line (not shown in Fig. ^ (b,d)) is only 
slightly affected by H. 

V. DISCUSSION 

We expect that our calculations are directly relevant 
to the quasi-2D layered materials La2-22:Sri+2a;Mn207. 
The existence of a strong superexchange interaction in 
this system is suggested by 

(i) relatively high values of Neel temperatures observed 
at the X = 1 endpointLj, 

(ii) intra-layer antiferromagnetic correlations present 
near Tc (Ref. |), 

(iii) interlayer (within the same bilayer) canting found 
at low temperaturegj (see also Ref. gj) and interlayer 
canted correlationsQ present above Tc- 

The latter point is associated with the structure of the 
quasi-2D manganates, which was discussed in the Intro- 
duction. 

The verification of the undersaturated behaviour ot 
low T remains an open question. It is not clear whether 
the materials La2-2i:Sri+2xMn207, with x — 0.4, lie 



within the region where the system exhibits undersatu- 
rated ferromagnetic behaviour at low T, or outside of this 
region (in the latter case, we still expect thermal fluctua- 
tions to be stronger than in a Heisenberg magnet, due to 
the presense of superexchange). Sqtoe measurements of 
the absolute value of magaetizatioiO in a; = 0.4 samples 
indicate undersaturationE3c3, while others do notcZl. 

We suggest that magnetic properties of the samples ex- 
hibiting undersaturation should be studied in the high- 
field, low-temperature regime of T ~ _ff . Our results (see 
Eqn. (|2^)) indicate that the intra-layer canted spin or- 
dering should be stabilized in this region. Another impor- 
tant prediction of our theory is the unusual dependence 
of the effective ferromagnetic exchange constant on mag- 
netization and hence on temperature (Eqn. (p^). While 
we did not study spin waves in the undersaturated low- 
temperature phase, it is clear that in such a situation 
the usual relationship between the low-T value of spin 
stiffness Dq and the Curie temperature {Dq oc Tc) is no 
longer valid. This might help explain the_recent experi- 
mental findings in perovskite manganatea^. We propose 
that the magnetization dependence of the effective ex- 
change constant (available through spin wave measure- 
ments) should be studied in more detail in both 3D and 
2D systems. 

It should be noted that the presence of undersaturation 
in ferro- and antiferromagnetic phases may well signal 
that in reality the system favours more complicated {e.g. 
spin glass-like, cf. Ref. spin ordering, that cannot 
be addressed within a single-site mean field theory. The 
fact that there have been no observations of an intra- 
layer spin canting in the layered compounds at, T ^ H 
is consistent with our results. 
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APPENDIX A: DERIVATION OF EQNS. (|-|). 

We begin by re-writing the local perturbation (||) as 

V ^ -{b-t){a\ai- ala2) (Al) 
in terms of the fermion operators. 
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11. . 

01,2 = ~q^i"^0) ± iC(l,0) + C(0,1) + C(_i,o) + C(0,-l)j 



2V2 



(A2) 



which anti-commute with each other. Perturbations of 
this form can be treated exactly by, [following 1. M. Lif- 
shits' theory of local perturbationsEjEJ. Here we will use 
mainly the Green's functions (resolvent operators) ap- 
proach of Refs. p^ , p^ . We will, without loss of generality, 
consider the 2D case. 



Perturbation (Al) results in a change of the net free 
energy of the carriers, which can be evaluated as 



SF = dn^ -T Tr <^ In 



-In 



1 + cxp 



T 



1 + exp 



f 



<y5(e) (i^tot(e) - Vtot{e)) dt . 



(A3) 



Here, ip{t) is defined by Eqn. (^, vtot{e) is the total (for 
the entire system ) ca rrier density of states in the presence 
of perturbation ([A1|), and 



AN 



7r2 2t- 



7IC 
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(A4) 



(where IC{x) is the complete elliptic integral and N is 
the number of lattice sites) is the total density of states 
corresponding to the unperturbed virtual-crystal band 
Hamiltonian, 



(A5) 



with the spectrum, e{q) — —t{cosqx + cosq^). The per- 
turbed density of states, i'tot(e), may include (5-function 
peaks corresponding to the discrete levels which split off 
downwards from the bottom of the band. As we shall see 
below, Oflly one discrete level may appear in the present 
problemLJ, and after integratio n Ip}^ . parts the free energy 
change ( |A3| ) can be re-writtenfijO in the form of Eqn. 
(^ (where the first term on the r. h. s. accounts for 
the contribution of the continuous part of the spectrum). 
We note that Eqn. (||) is an example of a Krein trace 
formula. The quantity 



(A6) 



is called the spectral shift function because of its relation- 
ship to the perturbation-induced shifts of the energy lev- 
els in thcj-aase of a discrete (or discretized) unperturbed 
spectrumllj. It can be evaluated as 

C(e) = -ImTr{lnG(e - iO) - lnGo(e - iO)} , (A7) 

TT 

where the operators Go(e) = (e • 1 — Hvc)^^ and 



G(e) = {G^\e) - Vy' ^ {1 - Go(e)y)"' Goie) (A8) 

are the Green's functions for the unperturbed and per- 
turbed Hamiltonians, respectively, and 1 is the identity 
operator. Eqn. ( [A7|) yields d^/de = -lniTr{G(e - iO) - 
Go(e — iO)}/7r. In turn, 

Tr(G-Go) -Tr{(i-GoF)-iGoFGo} = 
= TY{{i~G„V)-^GlV} = 

= 4 Tr ln(i - GoV) = 4 lnDct(i - GqF), (A9) 
de de 

where we used the fact that the operators (1 -pJjSaF)"^ 
and GqV commute with each other. ThereforeEJO, 



e(^) 



1 



ArgDet{i-Go(e-iO)y} 



(AlO) 



Since, according to Eqn. ( |Al[ ), the perturbation V is 
nothing but the sum of two projection operators, it is 
convenient to evaluate the r. h. s. of Eqn. ( AlO ) in a 
basis which includes the state s |1 ), |2), annihilated by 
the operators ai_2 (see Eqn. (A2)). In this basis, the 
determinant reduces to that of a 2 x 2 matrix, and one 
obtains 



e = Arg {1 + (& - t)iln - I22) -ib- tf x 

TT 

X [Det(/y ) - ^2 Det(G,,)] + Tri [{h - t)(Gii - G22)+ 
+ {b - tf{Gi2hi + G21/12 - Gii/22 - G22/11)] } . (All) 

Here the quantities lij and Cy /7r with i,j = 1, 2 denote, 
respectively, the real and imaginary parts of the matrix 
elements (i|Go(e — iO)|j). Explicitly, we find 



Cn(e) = l.(.)(l-^)^ 
Ci2(e) = G2i(6) = iKf) (1-^ 

We then obtain@ 



Det(Gy) = 0, Det(/y) = --2 



(A12) 



2t 

V 

-2t e 



Gijir])dT] 



2t 
2t 



v{7f}drj 



etc., and finally, Eqn. ( All ) takes form of Ec^n. 
latter can be conveniently re-written as 



(A13) 
I). The 



7rcot{7rf (e)} — — - 



1 



/(e) i2 _ 62 



i>{ri)djj 



(A14) 



where the branch of arc cot should be selected in a way 
which respects both the continuity of ^(e)lij and the fact 
that ^(e) = for 6 = 
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We note that the r. h. s. of the Eqn. (A14) has the 
form F(e)/i/(e), diverging as i^(e) — > 0. Therefore its 
possible values below the bottom of the band are Too, 
corresponding either to ^(e) = —1 or to ^(e) = 0. The 
case of i^(e) = —1 corresponds to the values of e between 
the bottom of the band and the bound state when the 
latter is present, while for e smaller than all the eigen- 
values of TC (continuous and discrete alike) the spectral 
shift function vanishes, ^(e) = 0. The change between 
these two values, which corresponds to the bound state, 
can occur only at e = zg, where zg satisfies the equation 
F{z) = 0. The latter condition yields Eqn. (|). Al- 
though it appears rather intuitive, this consideration of 
the bound-state problem can be substantiated by a direct 
calculation along the lines of Rcf. 

Interestingly, in Ref. |l^ the notion of a finite trace of 
certain operators in an infinite-dimensional Hilbert space 
(see Eqn. (A3)) was essentially introduced for the first 
time. Mathematical studies of related issues were initi- 
ated by M. G. Kreint^, and since then the Krein trace 
formulae remain an active research topic of functional 
analysis. 



APPENDIX B: SATURATED PHASES AT T=0 

In this Appendix, we present some numerical and ana- 
lytical results related to the phases which saturate at low 
T (i. e., the phases with values of magnetization or sub- 
lattice magnetization approaching unity at T — *■ 0). This 
being the simplest application of Eqns. (|^-||), it provides 
insight into the meaning of these equations which are cru- 
cial for the present paper. The undersaturated phases at 
low T will be considered in the Appendix 

We will assume that in the ground state, all the pairs 
of neighbouring spins form the same angle, 27, and that 
each spin forms an angle 7 with the z axis (thus, 7 = 
corresponds to the ferromagnetic phase, and 7 > - to 
the two-sublattice canted state of De Gennes) . Then the 
energy of the system at T = can be written ag3 



Eq \ cos "f + 2 Jaf cos 2'j — H cos J, (Bl) 



Here, Eq is the energy of electrons for i = 1. Let us now 
consider a single-spin perturbation of the ground state 
corresponding to the change of the polar angle value of a 
spin Si from 7 to a. The energy difference between this 
configuration and the ground state is given by 

dF^''\a)=SFDEibia),t,0) + 

+4:Jaf (cos(q; + 7) — cos 27) — (cos a — cos 7) , (B2) 

where (5F£is(6(a), i, 0) is given by Eqns. (§-^) with t = 
cos 7 and b ~ cos{(q! +iy)/2}. By minimizing with 
respect to cos 7, we findQ that the ferromagnetic phase is 
stable at I i?o I +H > 8Jaf- In this case one can use 
Eqn. (H) to obtain the value of SF^'^^ for a ^ 1: 



where the effective exchange constant (cf. Section III) is 
given by jj^^X = Jaf —\Eo\/8. In the case of pure double 
exchange, Jaf = 0, and for sufficiently large carrier con- 
centration X ~ 0.1, the effective exchange approximation 
is in fact adequate even for large values of a. Numer- 
ically, the difference between Eqns. (B2) and (B3) at 
a = TT does not exceed 15-20 %. This relative difi'er- 
ence (which reflects the different physics of the double 
exchange and Heisenberg exchange) becomes more pro- 
nounced at large Jaf ~ |-Eg|/8 (see Fig. ||). At smaller 
concentrations, x < 0.27 in 2D, and at sufficiently large 
values of Jaf, we find 5F'-^'>{n) < (dotted line in Fig. 

This means that the energy of the system can be 
lowered by hipping a single spin, and the ferromagnetic 
state becomes metastable. 

For larger values of Jaf, corresponding to Jaf > 
{\Eo \ + H)/8, the canted state with 



cos 7 



Eq I +H 
8Jaf 



(B4) 



emerges at T ^ 0, H > (see Sect. [II regarding the 
latter condition). In this case, the energy of small single- 
spin perturbations is quadratic in \b — t\/t <^ 1, 



SF^°\a) 



E, 







l&HJlp 



'AF 

' v'^M^\^^^(^e)de\{a-^fian^-1. (B5) 
2 - e / 



COS7- 



Note that the effective exchange approximation, which is 
based on the first-order (in {b — t)/t) perturbation theory 
result (^, is inapplicable. 

The typical results for SF'-^^{a) in the canted state 
are shown in Fig. ^ (left panel), where the dashed line 
represents the contribution of the band (first) term on 
the r. h. s. of Eqn. (||); one can see that the bound state 
noticeably lowers energies of fluctuations with a « 2tt—j. 

The origins of instabilities of the canted state which 
appear in our single-site treatment are illustrated in the 
right panel of Fig. ^ where functions SF^^'' (a) at differ- 
ent band fillings x for Jaf = 0.06, H — 0.01 are plotted. 
We see that as one lowers the bandfiUing from x — OA 
to X = 0.25, SF^^'>{'K — 7) becomes negative, so that the 
total energy can be lowered by flipping a single spin of 
sublattice I to the direction antiparallel to that of sub- 
lattice II spins, and the canted state is metastable. As 
one further lowers concentration to x = 0.15, the sign 
of d^5F^'^\a)/da'^ at a = 7 changes, signalling the in- 
stability of the canted phase. Indeed, since in 2D the 
principal- value integral on the r. h. s. of Eqn. (B5) di- 
verges at e ^ —2, the prefactor in front of (a — 7)^ in 
Eqn. (B5) is negative at small x. At H ^ 0, this coef- 
ficien t ch anges sign at x « 0.215 (cf. Appendix and 
Eqn.dCl). 



(5F(°)(a) = -iui°f\ -H){1- cos a) , 



(B3) 
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APPENDIX C: THE UNDERSATURATED 
FERROMAGNETIC STATE AT LOW T 

In this Appendix, we present results on the breakdown 
of the effective exchange approximation and on the low- 
temperature stability of the undersaturated ferromagnet. 

At 7? = 0, the first term in the expansion of SFi (see 
Eqn. (0)) in powers of 5A'I — cosai — M, 



5Fi{M, T) = A(M, T)5M + B{M, T)(5Mf 



is proportional to the effective exchange constant, A — 
4J^y(M). If the temperature is not too low, this linear 
term (which generates the effective exchange approxima- 
tion) provides a qualitatively reasonable approximation 
for 5Fi (see Fig. |2]). Thus, the system behaves as a 
Heisenberg ferromagnet with an M-dependent exchange 
constant. As explained in Sect. pH for sufficiently large 
values of Jaf (see Eqn. (|lT|)), jfff decreases with de- 
creasing T so that |jJ^y(M(r))| - T. Within the effec- 
tive exchange approximation, Afg — M(T) oc T at T — ^ 0. 

The effective exchange approximation, however, breaks 
down at low T, when the second term on the r. h. s. of 
Eqn. (CI) becomes dominant. This situation (which is 
depicted in Fig. ^ is due to the fact that the coefficient 



B{M,T) 



4^3 



So 



e^Ue)y^{e)de\ (C2) 



does not vanish at M 
Eqn. (1), and 



Mq. In Eqn. (C2), t is given by 



-2 



(C3) 



It is easy to see that the ratio of the 1. h. s. of (C4) to 
the r. h. s. increases as the maximum of i'o(e) at e = 
becomes more pronounced. Let us consider the extreme 
case of a constant density of states, t'o(e) = 1/4, and 
calculate both sides of (C4) at x = 0.5. We find: 



E^\=\^ £^ eVo(e)^^o(e)rfe = I [l^^^ + I 



(CI) so that the condition ( |C^ ) indeed is valid. Numerical 



calculations show that in 2D, inequality (C4) holds for 
X > Xc ~ 0.215. We anticipate that the value of Xc 
in 3D is lower. We also expect that, similarly to the 
ferromagnetic or canted state at T = (see Appendix 
^), the undersaturated FM state at low T may become 
metastable at values of x slightly above Xc- 

We note that the similar stability conditions for the 
antiferromagnetic and ca nte d (at small H) phases also 
take the form of Eqn. (|C4|). These should be distin- 
guished from the weaker thermodynamic stability condi- 
tion, d^/dx > 0, mentioned in Section [IV} The latter 
condition (in the antiferromagnetic, canted, and under- 
saturated ferromagnetic phases at iJ, T — > 0) can be re- 
written as 



1 



|£^ol 



8 Jaf [i^oifJ-o] 



>0, 



(C6) 



and in 2D holds at a; > 0.165. In writing Eqn. (C6), we 
assumed that M Mq at T ^ 0; note that this may be 
incorrect whenever the inequality (C4) is violated. 



Ho and i'o(e) are the chemical potential and the density 
of states in the unrenormalized {t = 1) band. 

At Mo - M{T) < VT, the hnear in 5M term in Eqn. 
( |cT|) can be omitted altogether. We then find that 5M = 
corresponds to an energy minimum if 



Eo 



> 



^- Jo{€)vo{^)de ■ 



(C4) 



In this case, the fluctuations of cos a at low T are conflned 
to the vicinity of Mo (see Fig. ^). 

When the inequality (C4) is not satisfied, the undersat- 
urated FM phase is expected to become unstable at low 
T. It is easy to see that in any dimensionality £> > 1, the 
inequality ( |C4| ) is violated at x — > 0. This follows from 
the fact that, when e approaches the bottom of the band. 



eJo(e) = 1 ^ / T^— di] > 1 . 



-D 



■q-e 



(C5) 



On the oth er hand, in 2D or in higher dimensions, the 
inequality (p3) is always satisfied for sufficiently large x. 
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FIG. 1. Single-spin fluctuation in the ferromagnetic phase. 
The bold arrow represents the average magnetization, and the 
dashed lines correspond to the hopping amplitude 6, which 
differs from the background hopping value t (solid lines). 



FIG. 2. Magnetization vs. temperature in the ferromag- 
netic phase at _ff = 0, a; = 0.4, and Jaf ~ 0.06. The solid, 
dashed, and dotted lines correspond to the 2D DE- superex- 
change magnet, effective exchange approximation, and usual 
Heisenberg ferromagnet, respectively. 



FIG. 3. The behaviour of the sublattice (solid line) and net 
(dashed line) magnetizations in the canted state ai H — 0.01, 
X — 0.4, and Jaf ~ 0.06, in comparison with the mag- 
netization of the ferromagnetic state (dotted line). The 
dashed-dotted line represents the results for the canting angle, 
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FIG. 4. Mean field results for the case of strong superex- 
change, Jaf = 0.08, x — 0.3. The solid and dashed-dotted 
lines represent the result for the sublattice magnetization m 
and canting angle 7 for H — 0.01. The dotted line corre- 
sponds to the sublattice magnetization niAFM of Neel AFM 
phase ai H = 0. 



FIG. 5. Phase diagrams of the DE-superexchange magnet 
for Jaf = 0.06 at H = (a) and H = 0.01 (b), and for 
Jaf = 0.05 at H = (c) and H = 0.01 (d), showing the 
ferro-, antiferro- (flop-phase at H > 0), paramagnetic, and 
canted phases (FM, AFM, PM, and CM, respectively). The 
dashed line in Fig. ^ (a,c) denotes the boundary of the ther- 
modynamically unstable [d^/dx < 0) region, and a; = a;i in 
Fig. ^ (c,d) separates undersaturated (x < x\) and saturated 
low-temperature regimes at H = 0. The behaviour of the 
system is symmetric with respect to quarter-filling, x = 0.5. 



FIG. 6. Single-spin perturbation energy, 5F^°\a), in the 
ferromagnetic state at T = in zero field, for x = 0.4 and 
Jaf = 0.04. Solid line represents the exact result (see Eqn. 
( |B2| )), while the dashed line corresponds to the effective ex- 
change approximation, Eqn. (B3). The dotted line represents 
the exact resuh for 5F^°\a) at x = 0.15, Jaf = 0.025. 



FIG. 7. The function 5^'°' (a) in the canted state at T = 0. 
The left panel corresponds to a; = 0.4, J^-ff = 0.08, H = 0.01. 
The dashed line represents the band contribution. In the right 
panel, 5F'^°\a) is plotted for J^// = 0.06 and H = 0.01 at 
X = 0.4 (solid line), x = 0.25 (dashed line) and x = 0.15 
(dotted line). 



FIG. 8. The energy cost, 5Fi, of a single-site fiuctuation 
in the ferromagnetic case at low T (see Appendix ^). The 
pronounced minimum of SQ at the average value of the po- 
lar angle a represents a sharp difference from the effective 
exchange approximation (dashed line) and causes the fluctu- 
ations of a (but not of the azimuthal angle, /3) to freeze out 
at low T. The plot corresponds to ii" = 0, a; = 0.4, T = 0.002, 
and Jaf = 0.06. 
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